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systems, namely (1) a classical system of two coupled oscillators, and (2) the 
Manko-Novikov metric (1992) which is a perturbation of the Kerr metric (a general 
'^ I relativistic system). We find their simple periodic orbits, their characteristics and 

their stability. Then we find their ordered and chaotic domains. As the energy 
goes beyond the escape energy, most chaotic orbits escape. In the first case we 
consider escapes to infinity, while in the second case we emphasize escapes to the 
central "bumpy" black hole. When the energy reaches its escape value a particular 

P3 , family of periodic orbits reaches an infinite period and then the family disappears 

(the orbit escapes). As this family approaches termination it undergoes an infinity 
of equal period and double period bifurcations at transitions from stability to 
instability and vice versa. The bifurcating families continue to exist beyond the 
escape energy. We study the forms of the phase space for various energies, and 
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^D ' increases abruptly as the energy goes beyond the escape energy. 
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1 Introduction 

The escapes from a dynamical system refer to a problem of basic interest for 
dynamics: the problem of chaotic scattering. This problem has attracted the at- 
tention of many authors in the previous decades (som e typical references are by 
[ ChurchiU et al. (1975); Petit a n d Henon^ (1986): Bleher et al. I Q988): Eckha r dtl 

T ■ 

( 



1988h : IJung and Schqlz l(ll988h:lH6ncml(ll989h:IContopoulos^^ and Kaufmann I 



19921 ) : iott and Tel I (|l993|) : lBenet et al. I (|l996lll998tn 



In particular we found ( Contopoulos and Efstathiou 1 12004|) the forms of the 



escape regions and the escape rates in the simple dynamical system 



H = ^{i +y+x+y)-xy=h (1) 



for various values of the energy h above the escape energy hesc = 0.125. We 
found also regions of orbits that never escape, or escape after a very long time. 
In particular the stable periodic orbits are surrounded by islands of stability that 
never escape. This problem is typical of "escapes to infinity". In fact, similar 
results were found in general polynomial Hamiltonians representing two perturbed 
harmonic oscillators. 

In the present paper we explore the escapes in the classical system Q and 
a very different problem of escapes, using a system from the domain of General 
Relativity. In particu lar the system we use belon gs to the so-called Manko-Novikov 
(MN) metric family. iManko and Novikovl ([1992') found an exact vacuum solution 
of Einstein's equations which describes a stationary, axisymmetric, and asymp- 
totically flat spacetime with arbitrary mass-multipole moments. The MN metric 
subclass that we use can be considered as a perturbation of the Kerr metric ( Kerr I 



subclass tnat we use can be consiaer eg as a p erturbation oi tne J\err metric i|J\err I 
Il963i) and it was introduced by Gair et al. I ([2008). The Kerr metric represents a 
rotating black hole of mass M and spin S. The MN perturbation is expressed by 
a parameter g, which measures how much the MN quadrupole moment Q departs 
from the Kerr quadrupole moment Qkctt = ^S^ jM (i.e. q = {Qkctt ^ Q)/M^)- 
While the Kerr metric describes an integrable system, the MN metric describes 
a non-integrable system that allows chaos. The line element of the MN metric 
expressed in the Weyl-Papapetrou cylindrical coordinates (p, (p, z) is of the form 



ds^ = - f{dt - Lodifif + r'^[e^^ {dp^ + dz^) + p^ dip^], (2) 



where /, oj, 7 are considered as functions of the prolate spheroidal coordinates 
V, w, while the coordinates p, z can be expressed as functions of v, w as well. Thus 



k^iv'^ — 1)(1 — w^), z = kvw (3) 



Periodic Orbits and Escapes in Dynamical Systems 



and 



„27 



e ^- 

2ke 

.27' 



B' 

2tpC 



4fc- 



A 

(t;2-l)(l-a2)2' 



^ = (t,2 _ 1)(1 + a bf - (1 - w^){b - ay, 
B=[{v+l) + {v- l)a bf + [(1 + w)a + (1 - 
C = (i.^ - 1)(1 + a b)[{b - a) - w{a + b)\ 

+ (1 - w^){b - a)[(l + a 6) + i;(l - a 6)], 
; «-P2 



.)6]^ 



In 



- w 
2 






+ 7TTTk(^3 ~ -P2) 



+ M-2+E 



V-W + {-lY^\v + w) 



£=0 



-aexp 



2 



-2/^ -1 + E 



^=0 



(i; - w)Pi 






Qexp 



2/?(l + E 






V «2 + «;2 



1, 



P. = P.(^). 



Here PeXC) is the Legendre polynomial of order I 



l)^ 



(4a) 

(4b) 

(4c) 

(4d) 
(4e) 

(4f) 
(4g) 



(4h) 
(4i) 

(4j) 

(4k) 
(41) 

(5) 



while the parameters k, a, /3 are related to the mass M, the spin 5*, and the 
quadrupole deviation q through the expressions 



-i+yi-x^ 



,k = M{^,P 






(6) 



while X is the dimensionless spin parameter x = S/M^. These formulae give the 
Kerr metric when q = 0. 

Contrary to the system ([T]), the escapes in the MN system refer not only to 
infinity, but al so to orbits t hat fall into the central compact object called "bumpy 
black hole" bv lGair et al. I (12008). The MN central compact object is not really a 
black hole, because its horizon is broken along the equator and regions of closed 
timelike curves exist outside the horizon (JGair et al. 1120081) . The MN problem is 
very different from the simple Hamiltonian ([T]!; however, the two problems have 
several common properties concerning the role of the periodic orbits and their 
relations to escapes. 
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In the present paper we study in detail the periodic orbits in the systems 
Q and ([2]) and their connection with the escape phenomena. We found that the 
periodic orbits in general do not lead to escapes, except in very special cases. 
Most of the families of periodic orbits generate an infinity of unstable, but non- 
escaping, periodic orbits. On the other hand, the asymptotic curves of the various 
unstable periodic orbits play an important role in generating chaos and most 
chaotic domains lead to escapes when the energy goes beyond the escape energy. 

The paper is organized as follows. In section[2]we calculate the periodic orbits of 
the system ([T]) . Then we find the asymptotic curves of the unstable periodic orbits 
and the corresponding chaos. In section [3] we find the escapes from the system 
([T]) and their statistics. In section |1] we find the periodic orbits in the system Q 
and the chaotic domains. In section [5] we find the corresponding escapes and their 
statistics. Finally, in section [6] we compare the two systems and draw our main 
conclusions. 



2 Periodic Orbits in the system (flj) 
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Fig. 1 (Tk) Simple periodic orbits (of period 1) for energy h = 0.12. 1^3) A surface of section 
{x,x) for h = 0.12 containing ordered and chaotic orbits. On the x = axis there are two 
stable periodic orbits (fe) and an unstable orbit (c) (at the intersection of the red and black 
curves), together with one stable (red) and one unstable (thick black) asymptotic manifolds. 
There are also two stable orbits (a) on the x = axis and triple and double period stable 
orbits that have bifurcated from the orbits (a) and (a') on the x=0 axis. 



For small values of the energy the system ([T]) has 4 periodic orbits of period 1. 
There are two straight line periodic orbits 



±V2 : 



(7) 



that are stable (orbits (a) and (a') in Fig. [T^). There is also a stable orbit (b), 
like an ellipse, described clockwise or counterclockwise, a simple unstable orbit 
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crossing the x — axis perpendicularly (orbit (c)), and the orbit y = (orbit (d)). 
In Fig. [TJj we see also 2 stable triple orbits 3 : 1 and 2 stable double orbits 2 : 1 
that have bifurcated from the orbits (a) and (a'). 

The orbits y = ±V2x undergo an infinity of transitions to instability and 
stability as h tends to hesc wh ile the period of these orbits tends to infinity 
( Contopoulos and Zikides 1119801) . The periodic orbits a and a terminate at h = 
hesc when the ir period is inf inite, in accordance with the Stromgren termina- 
tion principle ( Szebehelv Ill967il . For h > hesc the Curve of Zero Velocity (CZV) 
(a;^ + y^ —2 X y'^ = 2h) opens and the orbits a and a' escape to infinity. 

The families (a) and (a') at every transition to stability or instability generate 
by bifurcation new families of equal or double period. Families of higher order are 
also bifurcated from the families (a) and (a'). 

On a surface of section {x,x) with {y = and y > 0) the stable orbits are 
surrounded by islands of stability, while the unstable orbit is surrounded by a set 
of chaotic orbits (Fig. [1)3). In particular the orbits (a) and (a') appear on the axis 
a; = of Fig.[T]3. Near the unstable orbit (c) there is a chaotic region (Fig. [TJs). 

As h increases the chaotic region around the unstable orbit (c) increases, while 
various sets of higher order periodic orbits bifurcate from the straight line orbits 
{{a), {a')) [Fig.Wp)- 
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Fig. 2 Characteristics of some of the most simple families of periodic orbits. The family (a) 
and its bifurcations have i 7^ 



The characteristics of the various families of periodic orbits are shown in Fig. 
[2] These characteristics give the value of x when y = and y > 0. We mark in 
particular the orbits 3:1,2:1 and 1 : 1 bifurcating from the orbits (a) and (a) 
as the energy h increases. The stability of these families is given in Fig. [3l The 
orbits are stable if their Henon parameter a is between —1 and +1. 

The family y = \/2 x (a) is stable for < h < 0.111, then unstable in the 
interval 0.111 < h < 0.116, stable again for 0.116 < h < 0.1236, and has infinite 
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Fig. 3 SSh) The Henon stability parameter a of the families of Fig. [2] as a function of the 
energy h; ((Sji) Details of Fig. |3^. 



more transitions to instability and stability until h = hesc = 0.125. At the first 
transition to instability (h = 0.111) the value of a goes below a = — 1 and one 
stable family of double period bifurcates there. Its stability is represented by the 
red curve in Fig. [3l starting at a = 1. 

At the first transition from instability to stability along the family (a) {h = 
0.116) a goes above a = — 1 and there is a bifurcation of an unstable family 2 : 1 
(blue curve above a = 1). 

At the second transition to instability the value of a goes above a = \ (h = 
0.1236) and two different stable families of equal period bifurcate there (black 1:1). 
The orbits of these families deviate from the straight line y = \/2 x . At the second 
transition to stability there are bifurcations of two unstable families (magenta 1:1), 
and so on. 

Symmetric bifurcations are formed also around the family (a'). There are also 
higher order bifurcations from the families (a) and (a'), e.g. the family 3 : 1 (Figs. 
[2]and[3ji). 

The families bifurcating from the orbit (a) consist of orbits that reach the CZV 
at two points. These families do not terminate for h larger than h = 0.125, but 
extend all the way to h = oo. An example of an orbit of such a family (2 : 1) is 
shown in Fig. [3^ for h = 0.2. This family is bifurcated from the family (a) for 
h = 0.111 and it is stable in the interval 0.111 < h < 0.1252. M h = 0.1252 the 
family 2 : 1 becomes unstable and generates a family with double period (i.e. 4 
times the original period). For h = 0.137 the family 2 : 1 crosses again the. a = —1 
axis (Fig. [3^) and becomes again stable generating an unstable double period 
family (i.e. 4 times the original period). Finally, at h = 0.142 the family 2 : 1 
crosses the a = 1 axis and becomes unstable for larger h. At this crossing point 
it generates an equal period stable family (2 times the original period). This also 
becomes unstable for a slightly larger h, followed by a cascade of period doubling 
bifurcations, that generate an infinity of unstable families, that exist for arbitrarily 
large h. 
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Fig. 4 ||4ji) Some periodic orbits for h = 0.2. Four periodic orbits of period 1, (6), (c), and the 
Lyapunov orbits Oi, O2, and a periodic orbit of period 2, belonging to the family bifurcating 
from the family (a). ||4j5) An (unstable) periodic orbit (b) for h = 5.0. 



When h > hesc = 0.125 the CZVs open above and below the axis (y = 0) (Fig. 
|l^,b) and several orbits escape to infinity. However, the periodic orbits of Figs. 
|3^,b do not escape. E.g. the orbit 2 : 1 in Fig. S^ that bifurcated from the orbit 
(a) at h = 0.111 does not escape, although the orbit (a) has escaped. In Fig. [IJa 
we show a non-escaping orbit (b) for h = 5.0 when the openings are very large. 

At every opening of the CZV there is an unstable orbit (called Lyapunov orbit) 
crossing it (Fig. BK)- Any orb i t, crossing a Lyapunov orbit outwards, escapes to 
infinity ( Churchill et al. lll975tlContopoulos and Efstathiou 1 ,2004). The regions of 
direct escape are shown in red in Fig. [5] for h = 0.132. Orbits starting in these re- 
gions escape directly, without any intersections with the axis y = 0- The boundaries 
of the direct escape regions are defined by the asymptotic curves of the Lyapunov 
orbits. These a symptotic curves make infinite rotations around the escape regions 
( Contopoulos I1I99O I. In that paper we described in detail the asymptotic curves 
from the Lyapunov orbits. We found that these asymptotic curves include an in- 
finity of branches that start and end with infinite rotations around the escape 
regions. The stable asymptotic curves of any other unstable periodic orbit leading 
to escapes approach the same limiting asymptotic curves, i.e. the boundaries of 
the red regions as in Fig. [5] For relatively small values of the energy h, however, 
most of the chaotic orbits escape after one or more intersections with the axis 
y = (Fig.©. 

The asymptotic curves of the unstable orbit (c) surround the red escape regions. 
Examples of such asymptotic curves are shown in Figs. [T)d and [5] The unstable 
asymptotic curves start in opposite directions and the stable asymptotic curves 
start also in opposite directions. In Figs. [TJd and [5] we give in black the unstable 
asymptotic curve that starts towards the right downwards and surrounds clockwise 
the lower red region. As it comes again close to the point (c) it makes several 
oscillations, up and down close to the point (c). 
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Fig. 5 Surface of section for h = 0.132. Tiie orbits in tlic red regions escape directly to 
infinity without any intersection witii tiic axis y = 0. Blue dot corresponds to the orbit (c). 



The stable asymptotic curve (red) starts to the left and downwards, surround- 
ing counterclockwise the red escape region and intersects the unstable asymptotic 
curve at several homoclinic points. The oscillations of the stable and unstable 
asymptotic curves are much larger in Fig.[5](/i = 0.132) than in Fig.[T]D {h = 0.12). 
But the main difference between Figs.[T]D and [5] is that in the first case there are 
no escapes, while in the second case most of the chaotic orbits escape to infin- 
ity. The orbits that escape directly without any intersection with the axis y = 0, 
are marked in red. The stable asymptotic curve (red curve) of the unstable peri- 
odic orbit (c) (blue dot in Fig. [5]) makes an infinity of rotations around the lower 
red region. This asymptotic curve never intersects the red regions. This happens 
because orbits starting on the stable asymptotic curve have further intersections 
along the same asymptotic curve, until they reach asymptotically the periodic or- 
bit. Therefore they cannot escape immediately, hence they cannot be found inside 
a red region. 

On the other hand, the unstable asymptotic curve (black curve) intersects the 
red region after a number of oscillations below and above the stable asymptotic 
curve (Fig. [5]). The stable and unstable asymptotic curves intersect at an infinity 
of homoclinic points. A few points (1 — 5) are marked in Fig.[5l The correspond- 
ing homoclinic orbits approach asymptotically the periodic orbit (c) as f — > oo 
and t — > — oo. The unstable asymptotic curve (black) makes outer and inner loops 
above and below the stable asymptotic curve. The areas of these lobes are equal. 
The homoclinic intersections approach closer and closer the periodic orbit (c) , their 
distances decreasing proportionally to 1/A, where A is the larger (positive) eigen- 
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value of the orbit (c). Thus the lengths of the lobes become longer proportionally 
to A. The value of A increases as the energy increases. 

The inner lower lobes make some rotations, back and forth, around the lower 
red escape region. The inner loop starting downwards from the homoclinic point 
5 makes 1/2 rotation around the red region and returns closer than the point 5 to 
the periodic orbit (c). The next inner loop makes about 1.5 rotations around the 
red region and the subsequent loop intersects the red region (Fig.[5])- 

Then the orbits passing through points of the black curve inside the red region, 
escape to infinity without any further intersections with the surface of section. 
As a point along the black curve approaches the red region its next image on the 
surface of section makes an infinity of rotations around the upper red region of Fig. 
[5] approaching asymptotically the b oundary of the red region. Th i s phen omenon 
is described in detail in the paper of IContopoulos and Efstathiou I ( 20041) . 



For values of h up to about 0.32 there are some stable families, generated from 
the stable family (6), of orbits starting on the y = axis. This family consists 
of almost elliptical orbits that are described either clockwise or counterclockwise 
(Figs. [T^,b). The family (fe) is stable up to about h = 0.304 and then it becomes 
unstable at a period doubling bifurcation (a = —1). The double period families 
bifurcating from (6) start at a = +1 and become unstable at another period 
doubling bifurcation (4 times the original period) for about h = 0.32. Thus the 
interval between the first and the second period doubling bifurcations is Ah = 
0.016. This is followed by a cascade of period doubling bifurcations that produce 
an infinity of unstable families, existing for arbitrarily large h. It is known that the 
intervals between successive period doubling bifurcat ions in conserva tive systems 
decrease by a universal (asymptotic) factor S = ^) (JEckmann Ill98lh . Therefore, 
the total interval Ah that contains higher order bifurcation of the family (b) is 
approximately A = f^ = 0.018, and beyond h = 0.304 + 0.018 = 0.322 there is 

8.2 

no stable periodic orbit generated by bifurcations from the family (b). 

When the orbit (6) has just become unstable there are two islands of stability 
around it. For h = 0.305 (Fig.[S]) the asymptotic curves of the orbit (6) surround 
the two islands, but they do not extend to large distances. In fact, for this value 
of h there are invariant curves surrounding the orbit (fe) and the islands around 
it. However, for still larger h the asymptotic curves of (fe) extend very far (Fig. 
[7K) surrounding the whole available area on the surface of section (Fig.[7|D). The 
invariant curves surrounding the point (b) for h = 0.305, have been destroyed for 
h = 0.31 and the asymptotic curves come close to the limiting curve x +x = 2h. 

The same pattern of transition to instability is followed by most stable periodic 
orbits. Namely most stable orbits become unstable at a period doubling bifurcation 
followed by a cascade of infinite period doublings that lead to an infinity of unstable 
periodic orbits. An exception are the families (a) (a'), which have a non-unive rsal 
bifurcation ratio S = 9.22 and lead to escapes ( ContoDOulos and Zikides Ill980|). A 



theore tical explanation of this particular bifurcation ratio was provided bv lHeggie I 
(|l983t ). A r eview of previous th eoretical and numerical work on bifurcations is 
provided bv IContopoulos I ( 2002|) . 
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Fig. 6 The periodic orbit (fe) on the surface of section for h = 0.305, when the periodic 
orbit has just become unstable, and two stable orbits that bifurcated from it. The asymptotic 
curves from the orbit (b) surround both islands of stability, but they do not extend to larger 
distances. 
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Fig. 7 The asymptotic curves (blue=unstable and red=stable) of the orbit (b) for h = 0.31 
extend well beyond the areas around the two islands (Fig. [7^) and go all around the available 
area on the surface of section (Fig. [TJj) . 



3 Escapes from the system ([T]) 



As the energy h increases the proportion of orbits, that escape directly through 
the openings of the the CZV increases. In Fig. [5] (ft = 0.132) the upper red region 
represents orbits that escape through the upper opening of the CZV, and the lower 
red region represents orbits that escape through the lower opening of the CZV, 
without ever intersecting the axis t/ = 0. However, most of the chaotic orbits, 



Periodic Orbits and Escapes in Dynamical Systems 



11 



represented by scattered dots in Fig. [5l escape from the system after a small or 
large number of intersections with the y = axis. 
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Fig. 8 Escape regions for h = 0.15. 



As h increases to h = 0.15 the red regions increase in relative size (proportion 
of the total area). In Fig.[8]are shown the direct escape regions for h = 0.15 and the 
asymptotic curves (black, unstable and red, stable), from the periodic orbit (c), 
in the lower part of the figure. The regions in Fig. [8] are colored according to the 
number of intersections an orbit has with the axis y = before the orbit escapes. 
The stable asymptotic curve (red) again does not intersect the red region, but the 
(black) unstable asymptotic curve intersects the red region on the left after half a 
rotation around it. Then it intersects it again (after an inner and an outer loop, 
below and above the stable asymptotic curve) on the right side of the red region 
and so on. The orbits starting at points of the black curve inside the red region 
escape again without any further intersections with the axis y = 0. 

For larger h (h = 0.20) (Fig. ^) and h = 0.355 (Fig. [Qjs) the two regions of 
direct escapes increase in proportion to the total available space. 

In Fig. [9}3 we have marked one unstable asymptotic curve of the unstable 
periodic orbit (c) (black) , and one stable asymptotic curve of the unstable periodic 
orbit (b) on the left (red). The two curves are almost tangent at a point. We 
have checked that for h = 0.354 the two asymptotic curves do not intersect and 
for h = 0.356 they clearly intersect close to this point. Thus, for a value of h 
a little smaller than h = 0.355 the two curves are tangent. Then according to 
the Newhouse theorem (1973 Il983t ) close to the tangency point there is a stable 
periodic orbit, generated at a tangent bifurcation. 
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Fig. 9 S^) Escape regions for h = 0.25 eitlier directly (red), or after a number of intersections 
of the y = axis upwards (y > 0), as indicated in the color scale. The white region represents 
an island of stability. We mark also an unstable asymptotic curve (U) from the periodic orbit 
(c). l(9j5) The same, as in Fig. [9^, for h = 0.355, together with an initial arc of the stable 
asymptotic curve of the orbit (6) on the left (red). 



In Figs. Oi,b we see that the periodic orbits (6) and (c) are in a region where 
escapes occur after several intersections of the axis y = by the orbits. In fact, 
the periodic orbits never escape, and orbits close to them escape after a time 
that increases to infinity as the initial conditions approach the periodic orbits. In 
particular the stable asymptotic curves do not intersect the red regions of direct 
escapes. In fact, all the orbits starting on the stable asymptotic curves never escape. 

On the other hand, the unstable asymptotic curve U of (c) in Figs. [9^, b enters 
the red region at a rather small distance from the orbit (c), and the orbits starting 
inside the red region escape directly. However, this unstable asymptotic curve 
comes out of the red region several times, hence the corresponding orbits escape 
in general after some intersections with the axis y=0. 

The fact that there are chaotic orbits that escape after long ti mes is related to 

the phenom enon of stickiness along the unstable asymptotic curves ( Contopoulos and Harsoula I 
I2008ll2010[) . 

As h increases further new tangencies between the unstable and the stable 
asymptotic curves of the unstable orbits are formed, thus new stable periodic 
orbits appear. However, the islands around these periodic orbits are extremely 
small, since almost 100% of the orbits escape for large h. 

This situation is similar to the case where no escapes exist at al l, but tangencies 
betw een the asymptotic curves generate new stable periodic orbits ( Contopoulos et al. I 
Il994l' ). 



We proceed now to a statistical analysis of the escapes. This is done in four 
ways: 

(a) First we separate the orbits into ordered and chaotic + escaping. In Fig. llOh 
we give the proportion of the chaotic orbits N^i^g^/No (including the escaping 
orbits) as a function of the energy h. Namely we populate the whole circle inside 
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10a 



10b 




Fig. 10 I llOb ) The proportion of chaotic (including escaping) orbits N^-f^a/^O as a function 
of the energy h. mOb ) The proportion of the directly escaping orbits dNg/No (without any 
intersection with the axis y = 0) as a, function of the energy h. 
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The probability p„ 



of escapes as a function of n for various values of the energy h. 



the limiting curve x + x = 2hhy initial conditions of A'^o orbits (where A'^o = 
7r(250)^ ~ 200000) and we find the total numbers of chaotic + escaping orbits 
Ncha- Then the number of ordered orbits is A'^o — .^c/ia- When h is very small 
the proportion of ordered orbits is almost 100%. As h increases the proportion 
of chaotic orbits increases and for h = 0.125 this proportion is about 50% (Fig. 
llOh). When h goes beyond the escape perturbation h = 0.125 the proportion 
of chaotic orbits increases, but most chaotic orbits escape to infinity. Only 
a very small proportion of chaotic orbits near higher order unstable periodic 
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orbits inside islands of stability do not escape (like the orbits near the unstable 
periodic orbit (b) of Fig. [6}. 

The proportion of chaotic orbits (including the escaping orbits) can be given 
by the approximate formula (Fig. llOb ') 

N^ha / No =0.5[1+ tanh (30.0/1 -4.0)]. (8) 

For small h this proportion tends to zero, while for large h it tends to 1. 

(b) We calculate the proportion of direct escapes (Fig.llOb). dNo/No (where dNo is 
the number of orbits that escape without any intersection with the axis y = 0) 
as a function of the energy h. The direct escape regions are red in Figs. [5l 
[SI [9^,b. This proportion is small for h slightly larger than hesc = 0.125 and 
increases fast with h, reaching more than 90% for h = 0.7. 

(c) Then we study the escape time. The time is represented by the number of 
intersections n before escape. Most escapes take a long time for h slightly above 
hesc and take place faster as h increases. In Fig. Illh we give the logarithm 
of the proportion of escaping orbits d Nn/No where d Nn is the number of 
escapes after the nth intersection (upwards, y > 0) and before the (n + 1) 
intersection (i.e. there is no (n+ 1) intersection). We see that the escape rates 
are decreasing with n. The escape rates d Nn/No are large for small n and 
decrease considerably for larger n. For h > 0.35 the escapes take place quite 
fast. 

(d) Finally, we calculate the probability of escape as a function of time represented 
by the number of intersections n for various values of the energy. Namely the 
probability pn is the ratio d Nn/Nn where A''^ is the number of orbits that have 
not yet escaped before the nth intersection 

d Nn /„s 

This probability was calculated for the system ([1]) and other si mi lar sy stems in 
previous papers ( Contopoulos et al.lll993tlSiopis et al. Ill995ai bl. ll996h . It was 



found that pn tends to a constant value p, independent of the initial conditions 
of non-ordered orbits. For large times (large n) the probability of escape tends 
to zero if the energy is smaller than a critical value her, larger than the escape 
energy hesc- However, if h is larger than her, the probability p„ for large n, 
tends to a constant value p larger than zero. This quantity depends only on 
the energy h. In our present notation the limiting probability p is proportional 
to a power of the quantity {h — her), i.e. 

p Oc(/l-/lcr)", (10) 

where her is the critical value of h, and the exponent a is approximately a = 0.5. 
This value of a wa s found to be the same for different dynamical systems 
(jSiopiset al. lll995aHbh . 



In Fig. Illb we give the value of pn as a function of n for various values of h. 
When h = 0.15 the probability p„ tends to p = 0. On the other hand when 
h = 0.18,0.20,0.25,0.30 the probabilities tend to p = 0.006,0.018,0.026,0.034. 
These values can be approximated by the formula 

p = 0.11{h^hcr)°'^^. (11) 
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w here her = 0.175. This formula is consistent with the universal formula found 
bv lSiopis at al. I (|l995allbh . 



The statistics of escapes in various dynamical systems has been a subject of 
considerable interest in recent years. The main question is whether the various 
chaotic phenomena (like Poincare recurrences, corre lations and escapes) decay 
exponentially in time or according to power laws fe.g. lCristadoro and Ketzmerickl 
l2008t IVeneg eroles 2009^) . The power laws are supposed to be connected to the 
existence of stickiness in mixed dynamical systems that contain islands of sta- 
bility. However a detailed study of the effects of stickiness on the escapes should 
be the study of a future research. 



4 Periodic Orbits in the Manko-Novikov metric ((2) 

The Manko-Novikov (MN) spacetime depends on a parameter q that measures 
the quadrupole moment deviation of this metric from the Kerr metric with the 
same mass M and spin S. There have been prop osals and attempts to constrai n 
such deviations from observational data (see e.g. IPsaltis and Johannsen I ( 20101) : 



|Bambi (2011); Bambi and Barausse (2011a. b) and references therein). The final 
stages before the plunge ( escape) in a non-Kerr spacetime like the MN have certain 
astronomical interest fsee lBambi and Barausse I (|2011b|) ). From this point of view 
our work aims to exploit the non-integrability of the MN metric in order to give a 
detailed example of a such a final stage. 

The geodesic orbits of a test particle of mass /i are described as equations of 
motion of the following Lagrangian 

L = -^ ^. OtLv x'^x". (12) 

The MN metric has two integrals of motion ( Gair et al. l2008t iLukes-Gerakopoulos et al.l 



|2010[) . namely the energy (per unit mass) 

E = -^/l, = f{^^Lo^), (13) 

at 

and the z-component of the angular momentum (per unit mass) 

Lz = -^/fj- = fuj{i-uj ip) + f^^p'^ip, (14) 

wliere the dots mean derivatives with respect to the proper tim e. The Kerr m etric 
has one more integral of motion, the so-called Carter constant iCarterl ll968h . thus 
it is an integrable system. However, the MN model is non-integrable and it allows 
the appearance of chaos. 

The motion on a meridian axis {ip = cont) in the MN system satisfies the 
relation 

^{fP + i^) + V,ff{p,z)=Q, (15) 

where the effective potential Vf.ff{p,z) depends on q, E and Lz- Thus the motion 
takes place inside the CZV 



Veff - ie-- 



/-£"+ {^{Lz-i^E) 



0. (16) 
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Studies of the orbits in the MN metric we r e made by Gair et a l. (2008'l: lApostolatos et al. I 
( 2009h : lLukes-GerakoDoulos et al.n 20101 1 JContopoulos et al. i(|2011i') . In the present 
paper we study on detail the periodic orbits, chaos and escapes in this system. 

In the following we use the values q = 0.95, M = 1 and x = SjM^ = 0.9 for 
the spin. We find the orbits for a fixed value of the z-angular momentum Lz = 3, 
and various values of E. 



N 




Fig. 12 Three stable periodic orbits inside the outer CZV: mq, 2/3 and 2 = and two periodic 
orbits inside the inner CZV Uq (stable) and z = (unstable). There are also 5 CZVs that reach 
the axis p = 0. The line segment —0.43589 < z < 0.43589 on the axis p = represents the 
horizon of the bumpy black hole. The horizon is broken at the point p = z = 0. 



The form of the CZV in the case E = 0.95 is given in Fig. [121 It consists of two 
main closed curves and 5 more curves of small extent reaching the axis p = 0. The 
line segment —k < z < k {k = 0.43589 for M = 1, x = 0.9), along the axis p = 0, 
represents the horizon and any orbit reaching this axis escapes into the bumpy 
black hole. All orbits inside the 5 curves that are close to the horizon escape into 
the bumpy black hole. 

On the other hand, the orbits starting inside the main closed CZVs remain 
inside these curves for ever. In particular there are two simple periodic orbits 
inside the outer CZV, and two more simple periodic orbits inside the inner CZV. 
The periodic orbits mq (outer) and u'q (inner) intersect perpendicularly the p - 
axis while the other two periodic orbits are straight lines along the z = Q axis. 
Furthermore, we mark a stable resonant orbit 2/3. The orbits uq, u'q and the outer 
2 = orbit are stable, while the inner z = is unstable. 

If we take a surface of section (p, p) the intersections (z = 0, i > 0) of most 
orbits inside the outer CZV are along closed invariant curves around the point mq 
that represents the periodic orbit up (Fig. [T3b). However, there are also islands 
of stability around the stable (resonant) periodic orbits, and some chaos between 
these islands, around the unstable periodic orbits (Fig. 113b). 

On the other hand, most orbits in the inner region (inside the inner CZV) 
are chaotic (Fig. 113b ) although there is a large island of stability and three small 
islands of type 2/3. The chaotic character of these orbits is related to the instability 
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Fig. 13 Invariant curves and islands on a surface of section (p, p) {z = 0) for E = 0.95. I ll3b ) 
Islands and chaos in the inner region. I ll3b ) Invariant curves and islands in the outer region. 
The unstable periodic orbits are marked with big dots. 



I 




Fig. 14 Invariant curves around the periodic orbit Ug for E = 0.50. The boundary represents 
the orbit 2 = 0, which is stable in this case. 



of the orbit z = 0, which is represented by the boundary of the inner region on the 
surface of section. However, for a large range of values of E smaller than E = 0.67 
the orbit 2 = is stable and most orbits around u'q are ordered (e.g. for E = 0.5, 
Fig. [HI). 

If E increases beyond E = 0.9504 the inner and outer CZVs are joined into 
a common CZV (Fig. 115b . for E = 0.96). When the two curves join, an unstable 
periodic orbit is formed, which intersects perpendicularly the z = axis (point x in 
Fig.ll5b'l and exists for larger values of E. At the same time the two periodic orbits 
2 = join into one unstable orbit. The inner and outer regions are also joined on 
the surface of section (p, p) (Fig. 115b). The periodic orbit 2 = is represented now 



18 



George Contopoulos et al. 



15a 



15b 



8 


: 




/- — -^ 


. : 


6 


: 


y 


/ 


\ : 


4 


: 


czv/ 




\: 


2 





y 


Mo 

• 


[ 


-2 


: 


\ 




/: 


-4 


: 


\ 




/: 


-6 


: 


^ 


V 


/ : 


-8 


Till 




1 1 1 1 1 rrTTT-rrfTi 


1 1 1 1 1 1 1 i~ 



2 4 6 



8 10 12 14 16 18 

P 




-0.2 



-0.4 



8 10 12 14 16 18 
P 



Fig. 15 I llSt i) The two main CZV are joined into one for E = 0.96. Tlie point x represents 
an unstable orbit that is generated at the saddle point when the two curves join. I ll5b ) The 
corresponding surface of section (p, p). In this case the boundary, which represents the orbit 
2; = 0, is unstable. Chaos is dominant in the inner region and close to the boundary. 



by the common boundary of the two regions and close to it there is some chaos 
extending all the way around the orbit uq- Chaos is also dominant inside (i.e. on 
the left of) the orbit x fFig.lTSb'l. 
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Fig. 16 Characteristics of the periodic orbits {uq, x) and the boundaries of the main permis- 
sible regions (pis given as a function of the energy E) . 



Figure [16] gives the characteristics of the periodic orbits uq and x and the 
boundaries of the permitted motions for z = p = 0. The inner boundary is at 
an almost constant p ~ 0.72. The outer boundary increases considerably with 
increasing E and tends to infinity as _E — >• 1. When E > I the orbits escape to 
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infinity. When E < 0.9504 the permissible region splits into two and we have a 
middle outer boundary (MOB) (Fig. ll6|l which is the inner boundary of the outer 
region, and a middle inner boundary (MIB) which is the outer boundary of the 
inner region. 

As E decreases the outer region shrinks, and disappears for E ~ 0.9321. Then 
the orbit uq also disappears after becoming just one point. Then the orbit in 3 
dimensions is a circle on the equatorial axis around the center p = 0. For smaller 
E only the inner closed CZV persists. 

The position of the orbit u'q is close to the inner boundary, and cannot be 
distinguished from the inner boundary in the scale of Fig. 1161 In fact, the distance 
of this orbit from the inner boundary decreases, as E decreases, and goes to zero 
at about E = 0.28 fFig. llTb). For the same value of E the outer boundary of the 
inner region (MIB) reaches the inner boundary. Then the orbit uq in 3 dimensions 
is circular on the equatorial plane. 
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Fig. 17 1117b ) Periodic orbits u'q and x for E = 0.7. The family x was generated at the tran- 
sition of the orbit 2 = from stability to instability as E increases. Ill7b ) The characteristics 
of the families nj), x and the boundaries of the (inner) permissible region. 



As we have mentioned above the orbit 2 = of the inner region becomes stable 
when E decreases below E ~ 0.67. As E increases above this critical value, a 
stable family x that crosses perpendicularly the 2 = axis bifurcates from 2 = 
(Fig. llTb.b). As E increases further this family becomes unstable, generating by 
bifurcation higher order periodic orbits. E.g. for E = 0.95 one can see 3 small 
islands around this unstable family (on the right hand side of Fig. I13b l . 

When E decreases below E = 0.28 the inner CZV increases again (left side of 
Fig. 117b). and the orbit u'q deviates from the inner boundary. But, for E < 0.198 
the distance of u'q from the inner boundary decreases again fFig. llSp . although it 
does not tend to zero. 

The most strange evolution occurs when E becomes smaller than E = 0.19236. 
Then a couple of simple periodic orbits crossing perpendicularly the z = axis is 
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Fig. 18 Characteristics of the families ng,ng,nQ and of some resonant bifurcations from 
the family u'^ for relatively small values of the energy E. The dotted region indicates chaos. 
The dashed line gives the position of the first homoclinic point of the orbit Uq , and the gray 
region indicates escaping orbits, while the red region gives the initial conditions (for 2 = 0) of 
orbits that escape directly. 



formed out of nothing. One orbit (mq) is st able, and the other [ un 
lisp . This is called a "tangent bifurcation" ( Contopoulos Il2002r) . 



unstable (Fig. 



From the unstable point uq start two stable and two unstable asymptotic 
manifolds (Fig. ll9k l. The inner stable and unstable asymptotic manifolds surround 
the orbit up and intersect at an infinite number of homoclinic points. The first 
homoclinic point is on the axis p = 0, on the left of the orbit ug . The other couple 
of stable and unstable asymptotic curves approach the boundary and surround 
the orbit uq on the left of the figure. 

Near the unstable orbit Up and its homoclinic points there is some chaos (Figs. 
1181 and [T9b). As E decreases several periodic orbits bifurcate from u^ and near all 
the unstable orbits there are chaotic regions. These chaotic regions increase as E 
decreases. Two examples of bifurcating families, namely 2/3 and 1/3 are shown in 
Fig. [181 The orbit u^ remains stable from its generation at _E = 0.19236 down to 
E = 0.1912. Kt E = 0.1912 the orbit ug becomes unstable and a double period 1/2 
stable orbit bifurcates there towards smaller values of E. In Fig. 119b this orbit has 
receded from itg and two islands of stability are formed, one around each point. 
The orbits of the islands go alternatively from one island to the other. 



The orbit 



remains unstable for a small interval AE and then it becomes 



again stable, generating an unstable periodic orbit 1/2. 
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Fig. 19 ([T9ji) A surface of section (p, p) for E = 0.192 (just below the value of _E = 0.19236 
when the orbits «q , Uq are formed at a tangent bifurcation) . Invariant curves surround the 
orbits Mq , Uq and are also near the boundary. A little chaos appear near the unstable orbit u^ 
and its asymptotic curves. Ill9b ) A surface of section for E = 0.19045. It contains several islands 
around periodic orbits that have bifurcated from iig , like 1/2, 2/3 and 1/1. Several chaotic 
regions appear between the islands, but they are separated from each other by invariant curves 
around Mq . For this value of E the orbit Uq is unstable, while the boundary (representing the 
orbit 2 = 0) is stable. 



As E comes close lo E = 0.194045 the family itp becomes again unstable for a 
small interval AE. At the transition to instability it generates two different stable 
periodic orbits of equal period 1/1. In Fig. 119b the orbit uq is unstable and the 
two islands around it refer to different orbits. For a little smaller E the orbit u^ 
becomes again stable and generates two unstable periodic orbits 1/1. 

As E decreases further there is an infinity of transitions to instability and 
stability of orbits of double period 1/2 and of equal period 1/1. At the limit of 
this infinite sequence of transitions the orbit itg reaches the escape region (Fig. 
1181) and does not exist for smaller values of E. 



This phenomenon of infinite transitions to instability and stability along the 
same family was found also i n other dynamical systems, e.g. the system given by 
Eq. ([1]) of the present paper ( ContoDOulos and Zikides Ill980r) . 



The chaotic regions of Fig. 119b are separated from each other by invariant 
curves closing around uq . For a little smaller E most of these invariant curves are 
destroyed, as E decreases, and a large chaotic region is formed (Figs. [TSl and [20k ) 
for E = 0.19015. This large chaotic region contains several islands of stability, 
like 1/1, 1/2, 2/3 etc (Fig. 120b). that recede from Mq as E decreases (see the 
characteristics of the corresponding periodic orbits in Fig. I18p . 
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Fig. 20 Il20b ) A surface of section (p,p) for B = 0.19015. In this case most of the chaotic 
regions have joined. However, there are also islands, like 1/2, 2/3, 1/1. Il20b ) A surface of 
section for E = 0.1901 below the escape value of E = Eesc where escapes start to appear. The 
chaotic region of Fig. 120b is now a region of escapes, but we see also the same islands. 



5 Escapes from the Manko-Novikov system ([2]) 

In the MN system (including the Kerr system) there are two types of escapes of 
orbits: (a) Escapes to infinity when E > I, and (b) escapes to the central bumpy 
black hole for small values of E. 

The first type of escapes (escapes to infinity) is similar to the escapes in the 
Keplerian two-body problem (escapes along hyperbolic orbits). The permissible 
region increases in size as the energy increases. This size tends to infinity as the 
energy approaches the energy of the parabolic motions. 

On the other hand, the escapes to the bumpy black hole occur when the CZVs 
open inwards as the energy decreases and allow motions towards the bumpy black 
hole. 

In the MN case, that we study here, when E = 0.1901 the CZV of the inner 
region (the outer region does not exist) joins two CZVs, above and below it and 
orbits from the inner region are allowed to escape to the bumpy black hole (p = 0) . 
The aforementioned CZVs above and below the central region are the second and 
fourth leaf-like CZVs respectively (counting clockwise) in Fig. [12] from the set of 
five CZVs touching the broken horizon of the bumpy black hole [p = 0). These 
two CZVs expand outwards on both sides of the axis p = 0, while E decreases, 
until (for E « 0.1901) they reach the expanding central inner region (Fig. llSp . 
which till then contains only bounded orbits. After the connection, orbits of the 
inner region may escape downwards fFig.l21bl or upwards ('Fig.l21b) to the central 
bumpy black hole. 

In the MN case the escaping orbits are chaotic orbits that undergo chaotic 
scattering. In fact, if we compare Fig. 120b . {E = 0.19015) where no escapes are 
permitted, and Fig. 120b where we have many escapes we see a great similarity. 
That is the large chaotic domain around uq of Fig. I20fe i is transformed into an 
escape domain in Fig. 120b . We see several islands in Fig. 120b that are slightly 
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Fig. 21 Escapes in the MN metric Fig. 1121b ') downwards {E = 0.1901) and Fig. 1 121b ) upwards 
(E = 0.1899). 



changed from the islands of Fig. 120b . (the main change is that the islands 1/1 in 
Fig. 120b are surrounded by 3 secondary islands in Fig. I20bl . 

This similarity explains why the chaotic domain of Fig. 1181 (around uq ) is 
transformed into an escape domain for E < 0.1901. In particular the periodic 
orbit Uq of Fig. 120b reaches the point where the main CZV of the inner region 
joins the upper and lower CZVs and for smaller E this orbit escapes to the bumpy 
black hole, i.e. the periodic orbit does not exist any more. 
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Fig. 22 The period of the orbit u'q as a function of the energy E. 



This case is very similar to the cases y = ±^2a; (orbits (a) and (a')) of the 
system ([1]). There are two more similarities between the two cases. First is the 
sequence of infinite transitions to instability and stability of the orbits uq and the 
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orbits (a) and [a'). The second similarity is in the periods of the orbit Uq and of the 
orbits (a) and (a'), which both tend to infinity as the orbits approach the termina- 
tion point. The period of the orbits Uq is given in Fig. 1221 which increases consid- 
erably as E tends to _E = 0.1901 and presumably it tend s to infinity. The period of 
the o rbits y = ±^2x was also found to tend to infinity ( Contopoulos and Zikides I 
I198GI) . 

This behavior is consi stent with the S tromgren termination principle of the 
families of periodic orbits ( Szebehelv lllQGTTl . According to this principle the fami- 
lies either join other families, or terminate when their size, or their energy, or their 
period becomes infinite. 
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Fig. 23 1123b ) The proportion of the chaotic orbits (including the escape orbits), with initial 
conditions along the axis x, as a function of the energy E. Il23b ) The proportion dNo/No of 
direct escapes as a function of E. 



The proportion of chaotic -|- escaping orbits N^ji^/No increases as the energy 
decreases (Fig. 123b). We calculate the proportion N^i^^/Nq by taking initial condi- 
tions of orbits along the x-axis. This proportion tends to zero for relatively large 
E but it increases abruptly as E goes beyond Eesc (to smaller E). An approximate 
formula for this proportion is 



Ncha/No = 0.379[1 + tanh(304 - 1597S)] 



(17) 



and can be applied in the interval 0.188 < E < 0.192 (Fig. 25). The curve dTT]) is 
very similar to the corresponding curve of the Hamiltonian ([ij fFig.llOb.'l. However, 
for E < 0.186 the proportion of escaping orbits increases further and tends to 
100%. In Fig. 123b . two more abrupt increases can be seen. The first of these two 
increases takes place near E ~ 0.186, and around it one can use an approximate 
formula N^ha/No = 0.823 + 0.038 tanh(1399 - 7253S). The second increase takes 
place near E oi 0.183 and the corresponding approximate formula is N^ha/^o = 
0.921 -I- 0.065 tanh(163 - 892E) . This formula gives values near N^ha/No = 0.99 for 
E <0.180. 
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When E < Eesc most chaotic orbits escape to the bumpy black hole. However, 
only a proportion dNo/No escape directly from the system without any further in- 
tersection with the axis z = 0. This proportion increases as the energy E decreases 
(Fig. 123b V The curve dNo/No of Fig. 123b is very similar to the corresponding 
curve of Fig. llOb. In Fig. [18] we mark in red the initial conditions of the directly 
escaping orbits. Most of the other chaotic orbits escape after one or more intersec- 
tions beyond the original point, with the z = axis. But, as we see in Fig. 1181 for 
E < Eesc there are also orbits that do not escape. In particular the periodic orbit 
Ug and orbits close to it do not escape for much smaller values of E. Similarly the 
orbits that bifurcated from Uq, like the orbits 1/3, 1/2, 2/3, 1/1 (Fig. llSp exist 
also for E < Eesc- However, all these families become unstable for small E. In any 
case the existence of an infinity of families of unstable periodic orbits is a similar 
phenomenon as the one observed in the system Q. 



6 Conclusions 

We considered the periodic orbits and escapes in two quite different dynamical 
systems in order to find their common features. The first system is a system of 
two coupled oscillators, while the second system is a perturbation of the Kerr 
metric in General Relativity. 

In the system ([1]) escapes occur when the Curves of Zero Velocity open and 
many orbits escape to infinity. In the second case we have escapes to infinity, but 
we emphasized the escapes to the central bumpy black hole. In both cases when 
the energy changes and reaches the escape value, there is a periodic orbit that 
approaches the point where the CZV opens and its period tends to infinity. Thus, 
the family of periodic orbits terminates at the escape energy. Before the escape 
energy this family undergoes an infinity of transitions from stability to instability 
and vice- versa, and at every transition there is a bifurcation of an equal or double 
period family of periodic orbits, that does not escape as the energy varies. 

It seems that, in general, in systems with escapes there is a family of periodic 
orbits that terminates at infinite period when the CZVs open. 

Another common feature is that most chaotic orbits become escape orbits 
beyond the escape energy (beyond means larger energy in the first system, but 
smaller energy in the second system). However, for energies slightly beyond the 
escape energy most orbits take a rather long time to escape. As the energy increases 
(decreases) the proportion of fast escapes increases considerably. The proportion 
of chaotic (+ escaping) orbits increases abruptly as the energy goes beyond the 
escape energy. 

We studied the structure of the phase space and we distinguished the ordered 
and chaotic domains. The ordered domains surround the positions of the stable 
periodic orbits. Thus we found these periodic orbits and their bifurcations as the 
energy varies. On the other hand, the chaotic domains surround the asymptotic 
curves from the unstable periodic orbits. 

As the two systems under study are quite different we studied them in detail 
separately. 

In the first case (of two coupled oscillations) we found the characteristics of the 
periodic orbits and their stability. Then we found the structure of the phase space 
on a surface of section for various values of the energy. As the energy increases the 
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chaotic domains increase, while all stable periodic orbits become unstable. (How- 
ever, new stable periodic orbits appear for larger energies at tangent bifurcations.) 

The chaotic domains are covered by the asymptotic curves of a main unstable 
periodic orbit. As the energy goes beyond the escape energy most chaotic orbits 
escape, either directly, or after a small or a large number of intersections with 
the y=0 axis. Some asymptotic curves are split into successive pieces, each piece 
making infinite rotations around an escape domain. 

The statistics of chaos and escape are given by different graphs: (a) (b) the 
proportion of chaotic and escaping orbits and the proportion of the directly escap- 
ing orbits as functions of the energy, (c) (d) the proportion of the escaping orbits 
at the nth intersection with respect to the initial total number of orbits and with 
respect to the remaining orbits, as functions of the number of iterations (which 
represents time). 

In the system ^ (Manko-Novikov metric) we found various forms of the CZVs 
for various values of the energy, and the corresponding simple periodic orbits. 
We found also surfaces of section with mostly ordered orbits, for relatively large 
energies E (but E <1), or a mixture of ordered and chaotic orbits for smaller 
energies. We found the characteristics of the main families of periodic orbits and 
their bifurcations. 

We emphasized the appearance of a couple of periodic orbits (one stable, and 
one unstable) at a tangent bifurcation, and their evolution. As the energy decreases 
the stable family undergoes an infinity of bifurcations (at transitions from stability 
to instability and vice versa) . As the energy goes below a critical value the period 
of this orbit goes to infinity and then this orbit escapes to the bumpy black hole 
and this family of periodic orbits terminates. 

On the other hand, the bifurcating families exist and do not escape, but become 
unstable. We found the proportion of the chaotic regions as the energy decreases. 
Beyond the critical energy most of the chaotic region contains escape orbits. The 
proportion of the chaotic+escaping orbits increases abruptly below the escape 
energy and tends to 100% for smaller energies. 

It is remarkable that these quite different systems have very similar forms of the 
proportion of chaotic orbits fFigs. fTOk and 123b. ^ and of the proportion of directly 
escaping orbits (Figs.fTOb and 123b). The relations of the periodic orbits with the 
chaotic and escaping orbits are also quite similar. These similarities indicate that 
the properties of periodic, ordered, chaotic and escaping orbits that we studied in 
the present paper are very similar in generic dynamical systems. 
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